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and hydrogen-like ions in a superstrong magnetic field is studied. The finiteness of the nucleus 
size pushes the ground energy level up leading to a nontrivial dependence of the value of critical 



Qh nucleus charge on the external magnetic field. 
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I. INTRODUCTION 

At magnetic fields B > (37r/e 2 )B = 3iTml/e 3 w 6 • 10 16 GQthe Coulomb potential of the 
nucleus becomes screened due to large radiative corrections [T] (B = m 2 e /e = 4.4 ■ 10 13 G = 
4.4- 10 9 T). This leads in particular to the freezing of the ground state energy of the hydrogen 
atom at the value Eq = — 1.7 keV [U Ej. This statement is correct up to the values of the 
magnetic field at which the Landau radius an = 1/y/eB becomes close to the radius of the 
nucleus, a^ ps R. For hydrogen this happens at B « 10 5 B Q « 10 19 G, where the value 
of the proton charge radius R = 0.877 fm (see [3]) was used for numerical estimate. The 
approximation of a pointlike proton is not valid for B > 10 19 G, and in Sec. O and III we 
will find to what changes of atomic energies the finiteness of the proton size leads. 

With the growth of the nuclei charge Z the energy of the ground electron level diminishes 
and in the absence of magnetic field a nucleus with Z = 172 is critical: the ground level 
sinks into the lower continuum. As soon as the charge of the ion reaches Z cr = 172 two 
e + e~-pairs are produced from the vacuum. Electrons with opposite spins occupy the ground 
energy level while positrons are emitted to infinity [4j. In an external magnetic field the 
value of Z cr diminishes [HI E]. According to [5 J at B ~ 10 2 I?o uranium becomes critical, 
Z cr = 92, while at B = 10 4 i?o even comparatively light nuclei are critical, Z cr as 40. These 
results were obtained without taking the screening of the Coulomb potential into account. 
It was accounted for in [7J where it was found that because of screening larger magnetic 
fields are needed for a particular nucleus to become critical. Even more: according to [7] 
nuclei with Z < 50 do not reach criticality because of screening. The approximation of 
pointlike nuclei (which is valid when the Landau radius is considerably larger than the size 
of the nucleus) was used in [7J. A nucleus with Z w 56 becomes critical when the magnetic 
field is so large that the Landau radius equals the size of the nucleus. It was noted in [7J 
that the diminishing of the Coulomb potential due to the finite radius of the nucleus should 
push the value of the ground energy level up, and it was expected that this phenomenon will 



prevent ions with Z = 50, 51 from reaching criticality. In Sec. IV we present quantitative 
consideration on the influence of the finiteness of nucleus size on the dependence of Z cr on 
the value of the external magnetic field B. This dependence turns out to be not so simple. 
In particular, only nuclei with Z > 59 reach criticality in a superstrong magnetic fields. 



1 We use the system of units in which h — c= 1, a = e 2 = 1/137.03599 . 



When B further grows even nuclei with Z > 59 become noncritical (see Fig. pL 

In [U [2] and [7j the atomic nucleus was considered as an infinitely heavy source of the 
Coulomb field. Because of the finiteness of the mass of the nucleus its motion in the magnetic 
field should be taken into account and one should consider the two body (electron and 
nucleus) problem in the presence of a homogeneous magnetic field. This consideration was 
made in [8HTU] . According to the results obtained in [10] hydrogen atomic levels get increased 
by e\m\B/m p , where m is the projection of the relative angular momentum on the direction 
of the external magnetic field and m p is the proton mass. The corresponding formulae and 
numerical estimates are presented in Sec. |V} 

The finiteness of the nucleus mass also leads to a nonzero hyperfine interaction between 
the spins of the proton and the electron. Its importance in the case of a superstrong B was 



stressed in jTT] and we present our comments in Sec. VI 



II. ELECTRIC POTENTIAL OF THE NUCLEUS 

The finite size of the nucleus makes the Coulomb potential less singular at small distances, 
pushing up the electron levels. However the shape of potential at distances much smaller 
than the Bohr radius a# is not very important for the values of the electron energies. The 
effect is the largest for S levels, where the relative shift of energies goes like (R/clb) ~ 1CT 10 
where R is nucleus radius (in case of muonic atoms this shift is more important since it is 
enhanced by {m^/m e ) ~ 10 5 ). 

Strong magnetic fields make the Coulomb problem essentially one-dimensional. And in 
one space dimension a l/\z\ potential leads to a spectrum unbounded from below: the ground 
state energy equals minus infinity. The divergence of the potential at z — > is regulated by 
the Landau radius a#: \V{z)\ ~ e 2 /a#. It follows from this consideration that the behaviour 
of the potential at small distances determines the energy of the ground state. In this section 
we will find how accounting for the finite size of the nucleus modifies its electric potential. 

An analytic formula for the Coulomb potential of a pointlike charge along the z axis 
screened by a magnetic field was derived in [2j: 

$(0,z) = ^ (l - e^v 7 ^! + e-^ 1 ), (1) 



where z is the coordinate along the magnetic field, \x = x/Qml + (2e 3 B/n), and the charge 



is located at the point z — 0. The sum of the first two terms, ( 1 — e~' 2 ' V 6me I /|,2|, does 

not vary at distances z <C l/m e , and the smearing of the pointlike charge along the z axis 

within the domain with size R <C l/m e does not affect it. At the same time the last term, 

e~^ z \/\z\, becomes very sensitive to the charge distribution for B > l/(e 3 i? 2 ), jjR > 1. 

The potential in the plane transverse to the magnetic field (z = 0) was found in [7] for 

H > m e , B > m 2 /e 3 : 

6m 2 



$(p,0) = Mc-«»+-^pi, (2) 

where p is the coordinate in the transverse plane. The potential has a Yukawa behaviour 
both in the direction transverse to the magnetic field and along the z axis at distances 
p, z ^ Iq = - In £ — At these distances the potential is sensitive to the charge distribution. 



The important question is whether the long range part of the potential along the z axis, 
1 — e~' z ' V 6m e J /|^| ; is also affected by the smearing of the pointlike charge in the transverse 
plane. 

At z ^> l/m e the potential $(p, z) has the following behaviour [Tj: 

*0>,z) = -^= ^— rx ■ (3) 

' 2 +(l + ;0)p 2 



2 Z 



The main contribution from the terms ( 1 — e~' 2 ' v 6 ™ e ) /\z\ to the value of the electron energy 
comes from large distances, l/m e <C z ^C l/(e 2 m e ). To significantly change this contribution 
the potential ^ should noticeably differ from the Coulomb potential at distances z ~ 
l/(m e e): 

e3B R 2 > l ^ B > -^- w 4 ■ 10 10 B . (4) 

37rm 2 ~ (em e ) 2 ~ e 5 i? 2 

We are not going to consider here such strong fields; so this effect is neglected in what 
follows. 

Thus, there are two parts in the potential: the first one, ( 1 — e - ' 2 ' V 6 ™ e I /\z\, does not 
depend on the charge distribution inside the nucleus and the second one, originating from 
e -/w j r (where r = yfp 2 + z 2 ), is determined by the charge distribution. 

Another issue is the modification of the proton shape in a superstrong magnetic fields 
which would lead to the variation of atomic energies. As soon as the Landau radius of 
the electron becomes close to the proton radius the same happens with the Landau radius 



of the proton. When the magnetic field further grows one could expect that the size of 
the proton in the direction perpendicular to the magnetic field shrinks. But the proton is 
not an elementary particle and for R ^> a# the valence and sea quarks will be oscillating 
inside domains with size of the order of a# which are not necessarily situated at p = 0. The 
distribution of these rotating quarks inside the nucleus will be defined by strong interactions, 
so there is no reason to think that the nucleus in the transverse plane is squeezed down to 
the size of clhQ What is really happening to the proton shape in such a strong magnetic 
field is a subject for a separate study, while here we will neglect the possible shrinking of 
the nucleus in the magnetic field. 

With the account of the specific features discussed above the potential along the z axis 
at p ?S Ojj has the following form: 



1 - e - r v 6m l + h(R)e-» r , r>R, 
|(l- e- R v 6m - + h{r)e-» R , r < R, 



where h(r) is determined by the charge distribution inside the nucleus. 

Since the charge distribution inside the nucleus in such strong magnetic field is not known, 
we will consider three following cases: 

1. "Simple cut" - the potential outside the nucleus is equal to that of a pointlike charge 
and inside the nucleus it is constant and equals the value of the potential at the surface. 
This case corresponds to h(r) = 1: 



Ml- e~ r v 6m = + e-^ , r>R, 
$«(p^) = <| r V / ( 6 ) 

|(l- e~ R y 6m ' + e~f R , r < R. 



2. Homogeneously charged sphere — taking into account that the last term (p = eh(r)e ** /R 
in the expression (p3J) satisfies the Yukawa equation A(p — p?ip = — ATiQ{r) we obtain: 

■ (l _ e -rv^| + e -ixr . 1 ( ixR _ e -l*R\\ r >R, 

& 2 \p,z) = \ r \ ,— 2 " R[ V " (7) 

I ( 1 - e~ R y 6mi + e-» R ■ ^- (e^ r - e"^) , r < R. 



R \ ^ ' ^ 2/ir 



2 Nuclear core should prevent heavy ions from shrinking in the direction transverse to the magnetic field. 



3. Homogeneously charged ball — the potential can be easily found from the formula for 



a charged sphere: 

& 3 \p,z) = 



■ 1 1 _ e -V*"2 + e -» r ■ 2(^)3 (e^ R (pR - 1) + e'^^R + 1) ) ) , r > #, 



(8) 



III. THE FINITE SIZE OF THE PROTON AND THE HYDROGEN ATOMIC 
LEVELS 

The following equation for the hydrogen atomic energies on which the lowest Landau 
level (LLL) with m = splits was obtained in [12J by solving the Schrodinger equation (and 
it was checked in [7] that the relativistic corrections can be neglected as far as the binding 
energy is much smaller than the electron mass): 

21n^ + A + 21nA + 2^fl-^+4 7 + 21n2 = 2 [ dz [ l f°° (p)l „ d 2 p = /, (9) 
a B V A / J J J p l + z 1 

v r 

where the energies of the atomic states are determined by A, E = — (m e e 4 /2) A 2 , as = 
l/(m e e 2 ) is the Bohr radius, ip is the logarithmic derivative of the gamma function, and 



7 = 0.5772... is the Euler's constant. Roo(p) = e~ p / 4a » / ' J2iia 2 H is the wave function which 
corresponds to the transverse motion of the electron occupying the ground (n p = m — 0) 
Landau level. The dependence on the matching point z cancels in ^ for a H <C z <C 
l/(m e e 2 ). 

To take screening into account the factor 1/ 'v ' p 2 + z 2 in the right hand side of fel) should 
be substituted by: 



1 _», l (l _ p - V /p 2 +^ v /6^i +p - V /p 2 + 2 V6»nl+2 e 3i?/ 7 rA qqn 

Vp 2 + z 2 a/p 2 + z 2 \ ' ' 

which leads to the freezing of the values of atomic energies at B ^> m 2 /e 3 [TJ [2]. To take the 



finite proton size into account instead of (10) one should make the following substitution in 

@: 

' -+*(p,z)/e, (11) 



Vp 2 + z 2 



where $(p, z) is given by (J5J). For the right hand side of (J9J) we get: 

r VR 2 -z 2 
I = I 1 + I 2 + I 3 = 2jdz J 27rpdp \Roo(p)\ 2 ^ [l - e-V^f + h(y[? + z 2 )e 



-tiR 



+ 



R 



+ 2Jdz J 2irpdp\Rao(j))f 



VR 2 -z 2 

ZQ OO 



+ 2 [ dz f 2irpdp\Roo(p)f 



Vp 2 + z 



1 



V^T 






+ 



(12) 



R 

The sum I\ + I2 for a# ps i? is of order 1 and can be safely neglected. For an <C R we 
have: 



h + h « 2i2 • - 1 - e 



-R\/6m 2 



+ 



2 e -M« f R 



/" /i(s)dz = 2(l- e" 



R\/6m 2 



RV , R , --,- + (%""" ) -(13) 

We see that for three examples considered in Sec. |TT]the sum I1 + I2 rapidly diminishes when 
B grows, so we can safely neglect these termsjj 
For I3 we have: 



V R2+a ' 2 H 



hi 



dz 

z 

2(1 



1-e" 



'6m 2 , 



+ h(R)e 



/1Z 



R 2 + a 2 H 



Ex JR 2 + a 2 H JQm 2 e + h{R)Ex pJR 2 + a\ 



(14) 



where 



E x {x) = J 



t 



-dt. 



00 ( — X^rpr. 



n=l 



n • n\ 



E i( x )\ x ^i 



X 



- + 



1-2 1-2-3 



+ ... . 



;r-' 



(15) 



Substituting (14) in Q we get an equation which determines the values of A and the 
corresponding values of the energies on which the Landau level with n p = m = splits by 
the screened Coulomb potential: 



In— + ^ + \n\ + ^(l-\) -i 2- -1 In 2 
a B 2 V A/ 



In 



^0 



R 2 + a\ 



H 



Ex JR 2 + a 2 H J6m 2 + h{R)Ex pJR 2 + a%). 



(16) 



3 Let us note that the term (h)e tlR is not necessarily decreasing with the magnetic field because h(r) could 
correspond to any charge distribution inside the nucleus including a pointlike distribution. In case of a 



pointlikc distribution this term would prevent the ground energy level from going up. 



The dependence on the matching point zq cancels and finally we obtain the equation 
which determines the values of the freezing atomic energies with the account of the finite 
proton size: 

A 



- H-lnA + VU- v) H-2 T H-m2 = 

== - E x [^Jr 2 + a 2 H ^6m^j + h{R)E l [p^ R 2 + of) . (17) 



In 



In the limit B ^> 1/ (e 3 R 2 ) the right hand side of (17) does not depend on R and we 
obtain: 



1 

L - - 



A lim + 21nA lim + 2V'(l 



A 1 " 



+ 27 + 21n2 



In 



(18) 
(19) 



There are two ways to satisfy the equation ( 19 ) which has the large logarithm in the right 



hand side. The first one is to take a large A hm which will correspond to the ground energy 
level. The second one is to choose A lim close to the poles of ip{l — TiLr). The logarithmic 
derivative of the gamma function ip{x) has poles at x = 0,-1, —2 . . . which defines a series 
of A hm ~ 1/n, n = 1,2... This tower of A hm corresponds to the well known Balmer series 
of hydrogen atomic energies E^ m rj — (m e e 4 )/(2n 2 ). 

For the ground level from (19) we obtain A hm = 6.9 instead of the value for a pointlike 
proton A lim = 11.2 obtained in j2]. Let us note that for B ^> l/(e 3 R 2 ) all our approximations 
have a very good accuracy, so this result should have a good accuracy as well. To check it we 
solved the Schrodinger equation numerically. Since the adiabatic approximation is applicable 
(an <C a_B, 5> m 2 e z ) one has to solve the one dimensional Schrodinger equation with an 
effective potential V(z): 



dz 2 



2m e {E - V)x = 



(20) 



E 



m P e~ 



A 2 , V(z) = 



l H 



$(p, z) exp 



2a 2 Hi 



pdp. 



In order to numerically solve equation (20) an analytical formula for the averaged poten- 



tial energy V(z) is needed. For B ^> l/(eR 2 ), an <^ R, the potentials Q6J) — (JSJ) do not vary 
with p for p < a#, so the averaging over the transverse direction does not change them, 
i.e. V®(z) w -e$®(0,z), i = 1,2,3. For B < l/(ei? 2 ), a H > R, the modification of 



8 



the potential at distances r < R does not affect the electron motion while the averaging is 
very important since it removes 1/r singularity of $(p, z). So we need an analytical formula 
for the averaged potential energy for intermediate fields, B ~ l/(eR 2 ) at which an ~ R. 
However at distances r ~ R ~ an the screening does not occur (because 1/p « 10a#) so for 
B ~ l/(ei? 2 ) one should average the non-screened potential at these distances. 

Without taking screening into account the potential energy of the electron in the external 
electric potential of a homogeneously charged sphere has the following form: 



$ (0) (p,z) = I v^+- 



■, Vp 2 + z 2 >R, 






Vp 2 + z 2 < R. 

The formula for the corresponding averaged potential energy V^°'(z) looks like: 



(7(0) (. 



(i(i 



,(z 2 -R 2 )/2a 



I tt_ z 2 /2a 
a H V2 C 



erfc ( — 



V2 



2 — ,l^e z2 / 2a Heiic 



a H 



[*l 



«// 



v^ 






where erfc(x) is the complementary error function: 



(21) 



(22) 



erfc(x) 



erf(x) 



dy. 



(23) 



The extension of the formula (22) to the entire range of distances and magnetic fields for 



the potential energies V^'(z) and V^ 2 '(z) which correspond to the potentials <0>V-'(p,z) and 
$( 2 ^(p, z) is given by (i — 1, 2): 



V {l \z) 



2 c(0 / 1 



R 



1 



e (^-/? 2 )/2a 2 H 



-e 2 S[ 
-e^^/fe^Wcf--"- 



pE e z 2 /2a 2 H 
a H V 2 



«// 



>Off 



v/2 



51 



(1) 



l-e 

S{ 2) 



R\/6m 2 



+ e 



-fj.R 



, 5. 



(i) 

2 



erfc (^))' 1*1 <*> 



e -klv 6rn e + er^\ z \ 



I _ e " R V 6m e _|_ e "^ R 



2pkl 



,H z l 



-n\z\ 



s^ 



e -\z\y/6rn 2 _|_ g-/x|«| 



IpR 



(e^ R - e-^ R 



(24) 

(25) 
(26) 

(27) 



The formula (24) has the correct behaviour both for B <C l/(e 3 i? 2 ), ajj/e 3> -R, (because 
the screening factors S[ ~ 1 and Sv; ~ 1 for \z\ < R) and for B 3> l/(ei? 2 ), a# <C i?, 
(because averaging does not change &°\p, z): V ( -°\z)\B^i/( e R 2 ) ~ —e&°'(0,z)). Since these 



ranges of magnetic fields overlap the formula (24) is correct for all B. 



In Fig. [I] the behaviour of A gr (which corresponds to the ground level) according to the 

1 ( "Simple cut" ) and according to the results of the 



analytical expression (17) for h(r) 



numerical solution of Eq. (20) with V^ 1 '^) and V^(z) are shown. We see how A gr is going 



down (or the ground level is going up) when the Landau radius a# becomes of the order of 
the proton radius (B ~ l/(eR 2 ) ~ 2 • 10 5 Bq). The raising stops and the energy freezes at 
B ~ l/(e 3 R 2 ) « 3 ■ 10 7 B . According to Fig. |I|yW(^) and V {2) (z) lead to practically the 
same dependence of A gr on B so we use the potential V^ l \z) in what follows, n 

To present a qualitative explanation of the phenomenon of the rising of the ground energy 
level let us put zq « clb in Eq. rt9j). In case of a pointlike charge the main contributions to / 
(the right hand side of rt9J)) for B 3> m 2 e /e 3 come from integrating over Off < \z\ < l//x and 
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pointlike nucleus 

analytical (V (1) (2)) 

numerical (F (1) (2)) 

numerical {V {2) (2)) 
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10 1 



FIG. 1. The dependence of A gr on the magnetic field; the ground state energy equals Eq = 
— (m e e 4 /2) (A gr ) . The dot-dashed (red) line corresponds to the pointlike nucleus, the dashed 



(green) line — to the analytical formula (17), the solid (blue) and the dashed line with two dots 



(purple) — to the numerical solutions of (20) with V^(z) and V^ 2 '(z) correspondingly. 



4 The potential $' 3 '(p, z) varies inside the nucleus even in the non-screened case, that is why we do not have 
an analytical formula for the averaged potential energy. Nevertheless we made an estimate for V^(z) 
and found the numerical results for A gr which appeared to be rather close to the results obtained with 
VM{z) andi/< 2 )(z). 
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l/m e < \z\ < zq ~ as where the potential has the Coulomb (1/|^|) behaviour: 

J leB<i/^ « ln ^2 + ln (™e4) ~ In ^ + In - 4 = In -. (28) 

When the magnetic field grows the Landau radius an approaches the proton radius R and the 



first logarithm in (28) should be substituted by ln (l/(Rfi) 2 ) since at \z\ < R the potential 
does not have a l/\z\ behaviour. When B further grows and e 3 B approaches 1/R 2 the first 



logarithm in (28) goes away: 

AeSB^l/B? « ln (™e a l) = In ~ A - (29) 

The decreasing of / corresponds to the diminishing of A and, therefore, the ground energy 
level goes up. 

IV. CRITICAL NUCLEUS CHARGE 

It is well know that the Dirac equation in a pointlike Coulomb potential is not self- 
consistent for an electric charge Z > 137. If the finite size of the nucleus is taken into 
account then the Dirac equation becomes self-consistent. For Z rj 172 the ground energy 
level reaches the lower continuum, e = —m e , and two electron-positron pairs are created 
from the vacuum p[]. The electrons occupy the ground energy level while the positrons 
become free (two pairs are created due to the spin degeneracy of the ground energy level). 
This is known as the critical nucleus charge phenomenon. 

In the presence of a magnetic field the value of the critical nucleus charge diminishes [3] . 
The atomic energies were found by solving the Dirac equation. In [5j it was transformed 
into a set of two one-dimensional differential equation of first order which can be done when 
the adiabatic approximation is valid (B S> Bq (Ze 2 ) ): 

g z -(e + m e -V)f = Q, 
f z + {e-m e -V)g = 0, 

where g z = dg/dz, f z = dh/dz; e is the energy eigenvalue of the Dirac equation; the 
bispinor ip e = (^) of the electron is decomposed into ip e = l g{z)cxp (l p 2 /4a 2^), Xe = 

( if(z)ex (°- 2 /4a 2 ) J- The averaged potential energy V(z) is defined in the same way as it 
was done in Sec. IIII1 

11 



The analytical formula which describes the dependence of the atomic energies e on B was 
derived in [3]: 



Z, ' In | j V m ^_ £ 1 + ar ctan ( J™*±*\ + argT ( ^S + iZe 2 



-argr(l + 2zZe 2 )-^-(ln2 + 7 ) = ^ + mr , (31) 

Zi Zi 

where the argument of the gamma function is given by 

arg F(x + iy) = -yy + X) (r - arctan — — - j . (32) 

For the ground level at e > one should take n — 0, while for e < it should be changed 
to n = — 1. 



,,2 A 



Substituting £ = — m e into (31) the formula for the critical nucleus charge in a magnetic 
field was found in [5j: 

B nfr7 2 ^ 2 f , vr-2argr(l + 2iZ cr e 2 

According to this formula the critical nucleus charge diminishes with the magnetic field 
and for B w 10 2 i?o the uranium becomes critical (equations (31) and (33) are valid for 
B > max f (Ze 2 ) 2 B , B / (Ze 2 ^ 



2(Z^exp - 7+ «V . (33) 



To satisfy the matching condition used to derive the formula (31) the potential should 
be Coulomb at distances I > z Q , where z Q <C Ze 2 /(2m e ) is the matching point (see [5J for 
details). This condition is violated for the screened potential. Since we do not have an 
analytical formula for the energy levels which takes screening into account in the relativistic 
case, in [7] we solved the Dirac equation numerically. In order to do this we followed the 



paper [5J where the system (30) was transformed into one second order differential equation 

1 /9 

for g(z). By substituting g(z) = (e + m e — Vj x( z ) a Schrodinger-like equation for the 
function x( z ) was obtained in [5jjj 

d 2 X 



dz 



+ 2m e {E-U) X = , (34) 



E= £ ^i, U=^V-^V 2 + V " - + 3/8(y ' )2 



2m e m e 2m e 4m e (e + m e — V) m e (e + m e — V) 2 



5 This method of reduction of the Dirac equation to a Schrodinger-like equation was originally proposed by 
V.S. Popov to analyze the critical nucleus charge phenomenon qualitatively. 
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FIG. 2. The dependence of the ground energy level on the magnetic field for Z = 40. The dot- 
dashed (red) line corresponds to a pointlike potential without screening, the dashed (green) line 
- to a pointlike potential with screening, the solid (blue) line — to a potential which takes into 
account both screening and the finite size of the nucleus R = R^q ~ 5.1 fm. Cyan and purple lines 
(linestyles are shown in legend) corresponds to (hypothetical) larger nucleus radii, R = 3i?4o and 

it = 10it4Q. 



In [7] we used the equation (34) for numerical calculations. The numerical results for the 
pointlike nucleus were presented in [7] where the freezing of the ground energy level was 
obtained in the relativistic domain and the values of the critical magnetic fields were found 
for Z > 50. Nuclei with Z < 50 do not become critical due to screening. 

Our next step is to generalize the results of [7j in order to take the finite nucleus size 
into account. Let us note that in most cases the averaged potential energy V is smooth 
even when the potential $(p, z) has a cusp at the nucleus boundary. However due to our 
approximations in averaging V^'(z) done in the Sec. Ill the potential does have this cusp 



which is quite significant for B ^> l/(eR 2 ). It means that we have a ^-singularity in the 
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FIG. 3. The dependence of the ground energy level on magnetic field for Z = 40,59,60,90, 172. 
The correspondence between charge Z and color (linestyle) is shown in the legend. 

effective potential U and this should be taken into account in numerical calculations. 



Substituting V^ 1 '{z) for V(z) we checked that the equation (34) gives the values of the 
binding energy for hydrogen E = e ^™^ = — rn ^-\ 2 ^ which is very close to the ones obtained 
from the Schrodinger equation. 

Substituting the nucleus radius Rz = r v4 1//3 and V(z) = ZV^(z) into (34) we numer- 



ically calculate the value of the ground state energy of the hydrogen-like ion with charge 
Zq In Fig. [2] the dependence of the ground state energy on B is shown for Z = 40. We 
observe the rising of the ground energy level in the relativistic domain. The curves for a 
nucleus with Z = 40 and radii 3i? 40 and 10-R 40 are plotted to check how the energy depends 
on nucleus radius. We see that the rising starts at B ~ l/(eR 2 ) and stops at B ~ l/(e 3 R 2 ). 
The limiting energy does not depend on R. 

In Fig. [3] the dependence of the ground state energy on the magnetic field for Z = 
40,59,60,90,172 is shown. We see that ions with Z < 60 never become critical while 
a nucleus with Z = 60 is critical only within the small range of magnetic fields around 



In the numerical calculations we are using tq = l.lfm and A = 2.5Z. 
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formula (33), i?.=0 
with screening, R=0 
with screening, R^O 



critic ality 
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60 
Z„ 



FIG. 4. The values of magnetic fields at which the nuclei with charge Z cr becomes critical: a) 



without screening according to eq. (33), dot-dashed (red) line; b) numerical results with screening 



for pointlike nucleus, dashed (green) line; c) numerical results which take finite size of the nucleus 
into account, solid (blue) line. 



B Pd 10 4 i?o- For larger Z the range of magnetic fields in which ions are critical becomes 
wider. Ions become critical at B a little bit larger than the critical field for a pointlike 
nucleus. The rising of the ground energy level makes the ions noncritical for strong enough 
magnetic fields. Even an ion with Z = 172 becomes noncritical for B > 1.6 • 10 6 _Bn while 
it is critical for B = 0. To estimate the value of the nucleus charge at which the "final", 
or "second", freezing energy reaches the lower continuum we demanded that the nucleus 
should be critical for B = 10 8 So an d found that it is satisfied for Z > 210. It means that 
only ions with Z > 210 remains critical regardless of the value of the magnetic fields. 

In Fig. [4] the dependence of critical nucleus charge on B is shown. 
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V. THE HYDROGEN ATOMIC LEVELS IN A SUPERSTRONG B AND THE 
PROTON MOTION 

The two-body problem in the presence of a homogeneous magnetic field constant in time 
was analyzed in the papers [8HTU] where it was found that for an electrically neutral system 
the separation of the center-of-mass and relative coordinates can be carried out explicitly 
According to Eq. (12) from [10] the Hamiltonian which describes the relative motion of the 
electron and the proton in an external magnetic field B looks like: 



2M M V / 2m r 2 \ m P , 

(Bxr) + V(i 



e 2 f =. _\ 2 



8m r 

where the momentum K is an eigenvalue of the generalized momentum operator responsible 
for the motion of the center-of-mass of the atom, f= f e — f p and p — —id/dr are the relative 
coordinate and momentum of the electron and proton, M = m e + m p , m r = m e m p /M is 
the reduced mass, and V(r) describes the Coulomb interaction between the electron and the 
proton. The screening of the Coulomb interaction in a superstrong B should be taken into 
account as well [TJ |2] . 

Let us consider an atom at rest (K = 0) in a strong magnetic field B 3> m^e 3 in which 
the adiabatic approximation is applicable. The relative motion of the electron and proton in 
the direction of the magnetic field is determined by the potential V(z), while in the direction 
transverse to the magnetic field it is determined by B. 

Neglecting in a first approximation the Coulomb attraction for the atomic energy levels 



from (35) we obtain: 



,_, eB i \m\ + m + 1\ eB ( \m\ — m + 1 \ , , 

E ""» = ^ {"' + —2 ) + ^ ("' + —2 ) ' < 36 > 

In what follows we will be interested in the states of the hydrogen atom which originate 
from the lowest Landau level (LLL), for which n p = 0, m — 0, — 1, —2, . . . and the electron 



spin is antiparallel to B. For LLL the contribution of the first term in (36) summed with 
the energy of interaction of the electron spin with the external magnetic field B is zero while 
the second term gives a nonzero contribution. Thus the energies of the atomic states with 

different m are shifted by 

eB 

AE m = — [ml. (37) 

m p 
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Taking into account the motion along the z-axis governed by the (screened) Coulomb 
potential we obtain the following expression for the atomic energies: 

,_, m r e 4 , 9 eB\m\ ,_, eB\m\ ,„„ s 

£ = r —\ 2 + ^^ = E x + ^^, 38 

2 m p m p 

where in the first term m e can be safely substituted for the reduced mass m r since the 
numerical difference between m r and m e is very small. The values of A are determined by 
the following transcendental equation (see |2],Eq. (57)): 

ln( ^j— J =\ + 2\n\ + 2ip(l-\) + In 2 + 4 7 + ifj (1 + |m|) , (39) 

VI + ^/ V A7 

where H = B /(m? e e z ) is the magnetic field in units of atomic magnetic field, if) is the 
logarithmic derivation of the gamma function and 7 = 0.5772 ... is the Euler's constant. 



For each m — 0, — 1, —2, . . . the solution of (39) produces a tower of states; the tower with 
the lowest values of E\ corresponds to m = 0. The values of E\ in the limit B 3> 3nm 2 /e 3 
(or H 3> 37r/e 6 ) are shown in [2], Fig. 10. However the atomic energies of the states from 
different towers are shifted by the value AE m = eB\m\/m p and since AE m grows linearly 
with the magnetic field, for strong enough B this shift is big. 

Let us consider B = 2 • lO^m^e 3 ~ 4.7- 10 12 Gauss which was used in the calculation of the 



energies in the Table 1 of [10]. For m = for the ground level from (39) we obtain A[j = 4.3, 
E\{m = 0) = —255 eV. For m = —1 the ground level according to (39) corresponds to 
Aq" 1 = 3.8, E\(m = —1) = —193 eV which is 62 eV above E\(m = 0). The proton motion 
increase the atomic energies of the m — — 1 tower by eB/m p = 30 eV and the difference of 
energies with account of the finite proton mass becomes 92 eV (according to Table 1 of pH] 



it equals 93.3 eV). For stronger B the second term in (38) starts to dominate over the first 
one. 

VI. SPIN-SPIN INTERACTION IN HYDROGEN, HEAVY IONS AND POSITRO- 
NIUM IN A STRONG EXTERNAL MAGNETIC FIELD 

The interaction between the proton and electron spins leads to the hyperfine splitting of 
the hydrogen atomic levels. Being proportional to fi e fi p \ip(0)\ 2 ~ m e a 4 (m e /mj,) it is much 
smaller than the atomic energies, which are of the order of a 2 m e . In the case of a strong 
external magnetic field B ^$> m? e e 3 the spin-spin interaction considerably grows: 1^(0) Is ~ 
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|-0(O)| 2 f^ 8 -) = (^(O)! 2 (^^3)5 leading to a linear increase of the hyperfine splitting with 
magnetic fielq_] 

Ess-rriea 2 — — . (40) 

m p BQ 

It follows that at B ~ 10 5 i?o the spin-spin interaction energy becomes of the order of the 
freezing energy of the hydrogen ground level E$ m —1.7keV, and if a linear growth of Ess 
with B would take place for B > 10 5 Bo it would determine the value of the ground state 
atomic energy Eq. However just at B ~ 10 5 I?o the Landau radius an approaches the proton 
charge radius and a power formfactor suppression of Ess occurs preventing it from growing 
further. 

The energy of the spin-spin interaction in heavy hydrogenlike ions is enhanced by a factor 
Z which originates from the 1/as = Zm e a factor in the expression for |-?/>(0)| 2 . Since protons 
and neutrons from completely filled nuclei shells do not contribute to the magnetic moment 
of the nucleus, for B < 10 5 So the extra term in the energy is considerably smaller than the 
value of the electron mass and the consideration of nuclei criticality in strong B does not 
change substantially (One should also take into account that the formfactor suppression in 
heavy ions starts at smaller B). 

It would be very interesting to understand to which shift of the energy of the positronium 
ground state in a superstrong magnetic field the spin-spin interaction of the electron and 
positron leads. In the absence of external magnetic field, due to this interaction, the ground 
state of the parapositronium is lighter than the ground state of the ortopositronium 



E^S,) - E( l S ) = ^a 2 m e e 4 = ^e 8 m e . (41) 

The behaviour of the positronium energy levels in an external magnetic field has several 
specific features [H], and in fields B > e 8 ml/e the state with lower energy is a mixture of 
ortopositronium and parapositronium ground states in which the electron spin is oriented in 
the direction opposite to the magnetic field, while the spin of the positron is directed along 
B. Its energy shift due to the spin-spin interaction is of the order of: 

AE~^(0)\^ me e^, (42) 

B 



7 M. A. Andreichikov, B. O. Kerbikov, private communication 
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and for B > 10 4 i?o the ground state positronium energy could become lower than —2m e 
which should lead to the production of e + e~ pairs from the vacuum (see [15] as well). 

However the spin-spin interaction Hamiltonian is determined as a nonrelativistic expan- 
sion of the e + e~ scattering amplitude, the expansion parameter being p 2 /m 2 e [13] . In the 
case of a strong external magnetic field p 2 /m 2 ~ l/(a 2 H m 2 ) = B/B , and for B > B the 



correctness of the formulae (40) and (42) is doubtful. 

Let us remind that the anomalous electron magnetic moment leads to a linear growing 
with B of correction to the lowest Landau level energy. However this behaviour is valid only 
for B < B , while for B > B the strong linear dependence on B is replaced by a weak 
double logarithmic one ([IS], see also [2]). 

Concluding this section let us state that the behaviour of the spin-spin interaction in 
atoms and positronium at strong external magnetic fields B > Bq deserves further study. 



VII. CONCLUSIONS 

In a magnetic field B ^> m 2 /e 3 the potential of a pointlike charge becomes screened due to 
large radiative corrections. The screened potential has a Coulomb behaviour along the z axis 
at distances a H < \z\ < l/Ve 3 B and \z\ > l/m e and these regions define the ground state 
atomic energy. Distributing the pointlike charge within a domain of the size of the nuclear 
radius R leads to a less singular behaviour of the potential at distances \z\ < R. When e 3 B 
approaches 1/R 2 the ground state energy approaches the limiting value. For hydrogen the 
limiting value of A gr which defines the ground state energy Eq = — (m e e 4 /2)(A gr ) 2 is equal 
to 6.9 instead of 11.2 obtained in [2] for a pointlike charge, and Eq™ 1 = —0.65 keV instead 
of -1.7 keV. 

The same phenomenon of going up of the ground energy level is obtained numerically 
in the relativistic domain for hydrogen-like ions. It leads to a nontrivial dependence of the 
critical nucleus charge on the magnetic field: the nuclei with Z < 60 never become critical 
while the ions with 60 < Z < 210 are critical only within a finite range of magnetic fields. 
At Z = 210 a "second" freezing energy reaches the lower continuum and the nuclei with 
Z > 210 are critical at any B. 
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